N-Sulfonyl-1,2,3,4-tetrahydroisoquinoline Derivatives: Synthesis, Antimicrobial Evaluations, and Theoretical Insights.

Microbial contamination remains a significant economic challenge in the food industry, emphasizing the need for innovative antimicrobial solutions. In this study, we synthesized N-sulfonyl-1,2,3,4-tetrahydroisoquinolines (NSTHIQ) derivatives using an environmentally friendly Preyssler heteropolyacid catalyst, obtaining moderate to high yields (35-91 %) under mild conditions. Two derivatives (5 and 6) exhibited significant antifungal properties against various fungal species, including Aspergillus spp, Penicillium spp, and Botrytis cinerea. ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) analysis revealed the absence of hepatic toxicity in all compounds, making derivatives 2, 3, 4, and 5 potential candidates for further development. However, derivatives 6 and 7 exhibited immunotoxicity. In support of our experimental findings, reactivity indices were computed using Density Functional Theory principles, deriving valuable insights into the chemical properties of these derivatives. This study underscores the potential of NSTHIQ compounds as potent antifungal agents, coupled with the importance of employing environmentally friendly catalysts in drug discovery.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods.

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D2 dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Luminescent Lanthanide Metal Organic Frameworks as Chemosensing Platforms towards Agrochemicals and Cations.

Since the first studies of luminescent sensors based on metal organic frameworks (MOFs) about ten years ago, there has been an increased interest in the development of specific sensors towards cations, anions, explosives, small molecules, solvents, etc. However, the detection of toxic compounds related to agro-industry and nuclear activity is noticeably scarce or even non-existent. In this work, we report the synthesis and characterization of luminescent lanthanide-based MOFs (Ln-MOFs) with diverse crystalline architectures obtained by solvothermal methods. The luminescent properties of the lanthanides, and the hypersensitive transitions of Eu3+ (5D0→7F2) and Tb3+ (5D4→7F5) intrinsically found in the obtained MOFs in particular, were evaluated and employed as chemical sensors for agrochemical and cationic species. The limit of detection (LOD) of Tb-PSA MOFs (PSA = 2-phenylsuccinate) was 2.9 ppm for [UO22+] and 5.6 ppm for [Cu2+]. The variations of the 4f⁻4f spectral lines and the quenching/enhancement effects of the Ln-MOFs in the presence of the analytes were fully analyzed and discussed in terms of a combinatorial "host⁻guest" vibrational and "in-silico" interaction studies.

An inhibitory mechanism of action of coiled-coil peptides against type three secretion system from enteropathogenic Escherichia coli.

Human pathogenic gram-negative bacteria, such as enteropathogenic Escherichia coli (EPEC), rely on type III secretion systems (T3SS) to translocate virulence factors directly into host cells. The coiled-coil domains present in the structural proteins of T3SS are conformed by amphipathic alpha-helical structures that play an important role in the protein-protein interaction and are essential for the assembly of the translocation complex. To investigate the inhibitory capacity of these domains on the T3SS of EPEC, we synthesized peptides between 7 and 34 amino acids based on the coiled-coil domains of proteins that make up this secretion system. This analysis was performed through in vitro hemolysis assays by assessing the reduction of T3SS-dependent red blood cell lysis in the presence of the synthesized peptides. After confirming its inhibitory capacity, we performed molecular modeling assays using combined techniques, docking-molecular dynamic simulations, and quantum-mechanic calculations of the various peptide-protein complexes, to improve the affinity of the peptides to the target proteins selected from T3SS. These techniques allowed us to demonstrate that the peptides with greater inhibitory activity, directed against the coiled-coil domain of the C-terminal region of EspA, present favorable hydrophobic and hydrogen bond molecular interactions. Particularly, the hydrogen bond component is responsible for the stabilization of the peptide-protein complex. This study demonstrates that compounds targeting T3SS from pathogenic bacteria can indeed inhibit bacterial infection by presenting a higher specificity than broad-spectrum antibiotics. In turn, these peptides could be taken as initial structures to design and synthesize new compounds that mimic their inhibitory pharmacophoric pattern.

Rational Design of a Famotidine-Ibuprofen Coamorphous System: An Experimental and Theoretical Study.

Famotidine (FMT) and ibuprofen (IBU) were used as model drugs to obtain coamorphous systems, where the guanidine moiety of the antacid and the carboxylic group of the nonsteroidal anti-inflammatory drug could potentially participate in H-bonds leading to a given structural motif. The systems were prepared in 3:7, 1:1, and 7:3 FMT and IBU molar ratios, respectively. The latter two became amorphous after 180 min of comilling. FMT-IBU (1:1) exhibited a higher physical stability in assays at 4, 25, and 40 °C up to 60 days. Fourier transform infrared spectroscopy accounted for important modifications in the vibrational behavior of those functional groups, allowing us to ascribe the skill of 1:1 FMT-IBU for remaining amorphous to equimolar interactions between both components. Density functional theory calculations followed by quantum theory of atoms in molecules analysis were then conducted to support the presence of the expected FMT-IBU heterodimer with consequent formation of a R228 structural motif. The electron density (ρ) and its Laplacian (∇2ρ) values suggested a high strength of the specific intermolecular interactions. Molecular dynamics simulations to build an amorphous assembly, followed by radial distribution function analysis on the modeled phase were further employed. The results demonstrate that it is a feasible rational design of a coamorphous system, satisfactorily stabilized by molecular-level interactions leading to the expected motif.

Molecular design and synthesis of novel peptides from amphibians skin acting as inhibitors of cholinesterase enzymes.

Cholinesterases are a family of enzymes that catalyze the hydrolysis of neurotransmitter acetylcholine. There are two types of cholinesterases, acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), which differ in their distribution in the body. Currently, cholinesterase inhibitors (ChEI) represent the treatment of choice for Alzheimer's disease (AD). In this paper, we report the synthesis and inhibitory effect on both enzymes of four new peptides structurally related to P1-Hp-1971 (amphibian skin peptide found in our previous work. Sequence: TKPTLLGLPLGAGPAAGPGKR-NH2 ). The bioassay data and cytotoxicity test show that some of the compounds possess a significant AChE and BChE inhibition and no toxic effect. The present work demonstrates that diminution of the size of the original peptide could potentially result in new compounds with significant cholinesterase inhibition activity, although it appears that there is an optimal size for the sequence. We also conducted an exhaustive molecular modeling study to better understand the mechanism of action of these compounds by combining docking techniques with molecular dynamics simulations on BChE. This is the first report about amphibian peptides and the second one of natural peptides with ChE inhibitory activity. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

New small-size peptides modulators of the exosite of BACE1 obtained from a structure-based design.

We report here two new small-size peptides acting as modulators of the ß-site APP cleaving enzyme 1 (BACE1) exosite. Ac-YPYFDPL-NH2 and Ac-YPYDIPL-NH2 displayed a moderate but significant inhibitory effect on BACE1. These peptides were obtained from a molecular modeling study. By combining MD simulations with ab initio and DFT calculations, a simple and generally applicable procedure to evaluate the binding energies of small-size peptides interacting with the exosite of the BACE1 is reported here. The structural aspects obtained for the different complexes were analyzed providing a clear picture about the binding interactions of these peptides. These interactions have been investigated within the framework of the density functional theory and the quantum theory of atoms in molecules using a reduced model. Although the approach used here was traditionally applied to the study of noncovalent interactions in small molecules complexes in gas phase, we show, through in this work, that this methodology is also a very powerful tool for the study of biomolecular complexes, providing a very detailed description of the binding event of peptides modulators at the exosite of BACE1.

Looking for the interactions between omeprazole and amoxicillin in a disordered phase. An experimental and theoretical study.

In this paper, co-grinding mixtures of omeprazole-amoxicillin trihydrate (CGM samples) and omeprazole-anhydrous amoxicillin (CGMa samples) at 3:7, 1:1 and 7:3 molar ratios, respectively, were studied with the aim of obtaining a co-amorphous system and determining the potential intermolecular interactions. These systems were fully characterized by differential scanning calorimetry (DSC), FT-infrared spectroscopy (FTIR), X-ray powder diffraction (PXRD), scanning electron microscopy (SEM) and solid state Nuclear Magnetic Resonance (ssNMR). The co-grinding process was not useful to get a co-amorphous system but it led to obtaining the 1:1 CGMa disordered phase. Moreover, in this system both FTIR and ssNMR analysis strongly suggest intermolecular interactions between the sulfoxide group of omeprazole and the primary amine of amoxicillin anhydrous. The solubility measurements were performed in simulated gastric fluid (SGF) to prove the effect of the co-grinding process. Complementarily, we carried out density functional theory calculations (DFT) followed by quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analyses in order to shed some light on the principles that guide the possible formation of heterodimers at the molecular level, which are supported by spectroscopic experimental findings.

Structural and thermodynamic characteristics of the exosite binding pocket on the human BACE1: a molecular modeling approach.

We report a molecular modeling study aimed to locate and provide the full structural characteristics of the exosite binding site of the BACE1. A three-step procedure was followed. In the first stage, we performed blind docking studies on the whole target surface. In a second stage, the mode of binding was further refined by molecular dynamics (MD) simulation. Finally, binding free energy calculations, through the MM-PBSA protocol, were carried out to gain insight into the stability and thermodynamics of the inhibitor located at the selected binding pockets. Twelve binding pockets were identified on the surface of BACE1 by blind docking studies. The calculations of binding free energies for the 12 complexes show that van der Waals interactions dominate the mode of binding of these complexes. The best ranked complex shows that residues Glu255-Pro258, Phe261, Gly264-Ala272, Asp311-Ala313, Ser315, and Asp317-Tyr320 are located within 6 Å from the INH located at the exosite. The hydrogen bonds formed between the INH peptide, residues Tyr1, Tyr3, and Leu7 with the BACE1 residues Leu267, Cys269, Trp270, Asp311, and Asp 317 can strengthen the binding of the BACE1−INH complex.

Performance of density functional models to reproduce observed (13)C(alpha) chemical shifts of proteins in solution.

The purpose of this work is to test several density functional models (namely, OPBE, O3LYP, OPW91, BPW91, OB98, BPBE, B971, OLYP, PBE1PBE, and B3LYP) to determine their accuracy and speed for computing (13)C(alpha) chemical shifts in proteins. The test is applied to 10 NMR-derived conformations of the 76-residue alpha/beta protein ubiquitin (protein data bank id 1D3Z). With each functional, the (13)C(alpha) shielding was computed for 760 amino acid residues by using a combination of approaches that includes, but is not limited to, treating each amino acid X in the sequence as a terminally blocked tripeptide with the sequence Ac-GXG-NMe in the conformation of the regularized experimental protein structure. As computation of the (13)C(alpha) chemical shifts, not their shielding, is the main goal of this work, a computation of the (13)C(alpha) shielding of the reference, namely, tetramethylsilane, is investigated here and an effective and a computed tetramethylsilane shielding value for each of the functionals is provided. Despite observed small differences among all functionals tested, the results indicate that four of them, namely, OPBE, OPW91, OB98, and OLYP, provide the most accurate functionals with which to reproduce observed (13)C(alpha) chemical shifts of proteins in solution, and are among the faster ones. This study also provides evidence for the applicability of these functionals to proteins of any size or class, and for the validation of our previous results and conclusions, obtained from calculations with the slower B3LYP functional.

Predicting 13Calpha chemical shifts for validation of protein structures.

The (13)C(alpha) chemical shifts for 16,299 residues from 213 conformations of four proteins (experimentally determined by X-ray crystallography and Nuclear Magnetic Resonance methods) were computed by using a combination of approaches that includes, but is not limited to, the use of density functional theory. Initially, a validation test of this methodology was carried out by a detailed examination of the correlation between computed and observed (13)C(alpha) chemical shifts of 10,564 (of the 16,299) residues from 139 conformations of the human protein ubiquitin. The results of this validation test on ubiquitin show agreement with conclusions derived from computation of the chemical shifts at the ab initio Hartree-Fock level. Further, application of this methodology to 5,735 residues from 74 conformations of the three remaining proteins that differ in their number of amino acid residues, sequence and three-dimensional structure, together with a new scoring function, namely the conformationally averaged root-mean-square-deviation, enables us to: (a) offer a criterion for an accurate assessment of the quality of NMR-derived protein conformations; (b) examine whether X-ray or NMR-solved structures are better representations of the observed (13)C(alpha) chemical shifts in solution; (c) provide evidence indicating that the proposed methodology is more accurate than automated predictors for validation of protein structures; (d) shed light as to whether the agreement between computed and observed (13)C(alpha) chemical shifts is influenced by the identity of an amino acid residue or its location in the sequence; and (e) provide evidence confirming the presence of dynamics for proteins in solution, and hence showing that an ensemble of conformations is a better representation of the structure in solution than any single conformation.

New antitumoral acetogenin 'Guanacone type' derivatives: isolation and bioactivity. Molecular dynamics simulation of diacetyl-guanacone.

We describe herein the isolation and semisynthesis of four acetogenin derivatives (1-4) as well as their ability to inhibit the mitochondrial respiratory chain and several tumor cell lines. In addition, four nanoseconds (ns) of MD simulation of compound 4, in a fully hydrated POPC bilayer, is reported.

Position dependence of the 13C chemical shifts of alpha-helical model peptides. Fingerprint of the 20 naturally occurring amino acids.

The position dependence of the (13)C chemical shifts was investigated at the density functional level for alpha-helical model peptides represented by the sequence Ac-(Ala)(i)-X-(Ala)(j)-NH(2), where X represents any of the 20 naturally occurring amino acids, with 0 < or = i < or = 8 and i + j = 8. Adoption of the locally dense basis approach for the quantum chemical calculations enabled us to reduce the length of the chemical-shift calculations while maintaining good accuracy of the results. For the 20 naturally occurring amino acids in alpha-helices, there is (1) significant variability of the computed (13)C shielding as a function of both the guest residue (X) and the position along the sequence; for example, at the N terminus, the (13)C(alpha) and (13)C(beta) shieldings exhibit a uniform pattern of variation with respect to both the central or the C-terminal positions; (2) good agreement between computed and observed (13)C(alpha) and (13)C(beta) chemical shifts in the interior of the helix, with correlation coefficients of 0.98 and 0.99, respectively; for (13)C(alpha) chemical shifts, computed in the middle of the helix, only five residues, namely Asn, Asp, Ser, Thr, and Leu, exhibit chemical shifts beyond the observed standard deviation; and (3) better agreement for four of these residues (Asn, Asp, Ser, and Thr) only for the computed values of the (13)C(alpha) chemical shifts at the N terminus. The results indicate that (13)C(beta), but not (13)C(beta), chemical shifts are sensitive enough to reflect the propensities of some amino acids for specific positions within an alpha-helix, relative to the N and C termini of peptides and proteins.

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide.

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.

Polyproline II helix conformation in a proline-rich environment: a theoretical study.

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP(II) structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro)(3)-X-(Pro)(3)-Gly-Tyr-NH(2), with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro)(3)-Ala-Ala-(Pro)(3)-Gly-Tyr-NH(2) (PAAP), by using a theoretical approach that includes a solvent effect as well as cis <--> trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since (13)C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles (phi,psi) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged (13)C chemical shifts for the guest residues in the PXP peptide (X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.

Unblocked statistical-coil tetrapeptides in aqueous solution: quantum-chemical computation of the carbon-13 NMR chemical shifts.

We recently reported a theoretical characterization of representative ensembles of statistical-coil conformations for tetrapeptides with unblocked termini in aqueous solution, at pH 7. The results showed good agreement between the computed Boltzmann-averaged and experimentally-determined values for both the vicinal coupling constants (3)J(NHalpha) and the alpha-proton chemical shifts. Here, we carry out a cluster analysis of the ensembles of conformations generated in that study, and use them to compute the Boltzmann-averaged values of the quantum-chemical (13)C chemical shifts for different amino acids in the unblocked tetrapeptides GGXA (where X stands for Phe, Arg, His, Glu, Ile, Lys, Gln, Tyr, Leu, Thr, Ala, Gly and Val). The values of the (13)C chemical shifts in these thirteen amino acids (for which experimental data are available) were computed by using Density Functional Theory with a 6-311+G(2d,p) basis set. Good agreement is found in terms of both the correlation coefficient (R) and standard deviations of the difference between the computed Bolztmann-averaged and the NMR-determined values for the (13)C chemical shifts. These results suggest that it may be possible to build a reliable theoretically-derived database of chemical shifts for statistical-coil residues. The results of the current study contribute to our understanding of the relations between chemical shifts, dihedral angles and vicinal coupling constants, (3)J(NHalpha). In addition, they can shed light as to how the statistical-coil conformation is related to the conformational preference of more structured states, such as the alpha-helical conformation.

Unblocked statistical-coil tetrapeptides and pentapeptides in aqueous solution: a theoretical study.

NMR studies of the molecular conformations of peptides and proteins rely on a comparison of the relevant spectral parameters with the corresponding values for so-called statistical-coil polypeptides. For this reason, it is necessary to characterize the experimental ensemble of states populated by statistical-coil peptides. Such a characterization, however, has proven to be both difficult and sensitive to changes in many environmental parameters such as solvent composition, temperature, pH, as well as the neighboring amino acids in the sequence. As a consequence, a series of significant discrepancies has been reported for some experimentally observed parameters, such as chemical shifts, or vicinal coupling constants, (3)J(NHalpha), whose values appear to be incompatible with a statistical-coil ensemble. In this work, we report the results of a molecular mechanics study of a series of unblocked tetra- and pentapeptides under different pH conditions. These calculations were carried out with explicit consideration of both the coupling between the process of proton binding/release and conformation adopted by the molecule at a given pH and the contribution of the conformational entropy to the total free energy. Good agreement was found between the calculated and experimentally determined values of the vicinal coupling constant, (3)J(NHalpha), the alpha-proton chemical shift, and the (13)C(alpha) chemical shift. All the evidence accumulated in these theoretical calculations helps to rationalize some of the unsettled anomalies observed experimentally, and to provide an understanding of the effect of pH and amino acid sequence on the conformational preferences of statistical-coil peptides.